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Abstract 

We studied SU(2) flux distributions on four dimensional Euclidean lattices 
with one dimension very large. By choosing the time direction appropriately 
we can study physics in two cases: one is finite volume in the zero temperature 
limit, another is finite temperature in the intermediate to large volume limit. 
We found that for cases of j3 > (3 C there is no intrinsic string formation. Our 
lattices with (3 > j3 c belong to the intermediate volume region, and the string 
tension in this region is due to finite volume effects. In large volumes we found 
evidence for intrinsic string formation. 
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I. INTRODUCTION 



There are two complementary ways to obtain approximate solutions of strong coupling 
QCD. One is to take 3-space to be a finite volume torus and obtain semi-analytic solutions 
|l||JQ] of an effective Hamiltonian. The second is via lattice simulations. Confinement can 
be studied by each method in limited domains. In small physical volumes the fields are 
very rigid and the problem can be treated using a variational method applied to a small 
number of dynamical variables. In this domain the string tension is found to be zero. 
This result is consistent with asymptotic freedom since only short wavelength modes occur. 
At intermediate volumes there is a clear signal for string tension and further that it is a 
consequence of a tunneling amplitude between the vacuua that are degenerate for small 
volumes. In this domain, lattice methods are also accessible and are in good agreement for 
quantities such as glueball masses and string tension. Lattice calculations can take over to 
study larger volumes where semi-analytic methods become prohibitive. 

The existence of string tension in finite volumes does not imply confinement. Clearly if 
the the volume is not large enough to allow the fields to spread out, the finite box itself may 
be responsible for the linearly rising potential energy between quarks. Global studies have 
left open the question of the volume at which intrinsic confinement takes over ||. In this 
paper we look at a local quantity, the flux tube profile between static quarks as a function 
of physical volume in SU(2) lattice gauge theory in order to elucidate this question. 

The physical size of the box is characterized by a dimensionless variable 

z g = m(0+)L, (1) 

where m(0 + ) is the lowest glueball mass which is the energy gap or in terms of length it is 
the inverse of the correlation length. The length L is the linear size of the box. M. Luscher 
studied QCD in a small box, z g < 1, with periodic boundary conditions. He derived a 
low-energy effective Hamiltonian for SU(N) gauge theory in small volumes 0, i.e. z g < 1. 
Subsequently the lowest energy levels of SU(2) and SU(3) || gauge theories in small 
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volumes were computed by using this Hamiltonian. Van Baal and Koller then found that 
the crucial tunnelling between degenerate vacuua can be obtained by imposing appropriate 
nonperturbative boundary conditions on the Raleigh- Ritz trial wave functions 0]. They 
extended the calculation of the SU(2) glueball masses up to z g w 5.0. 

Berg and Billoire U carried out a thorough study of glueball masses, electric flux states 
and string tension for intermediate volumes (1 < z g < 5). They provided a detailed compar- 
ison between their numerical results and the analytic results of Van Baal and Koller. They 
chose lattice sizes N a x N£ x N c , with 

N a < N b < N c . (2) 

By identifying N c to be the time extent, one can simulate the zero temperature finite volume, 
N a x N£a 3 , field theory, with the temperature defined as || 

where Tb is called box temperature in ref. ||, and a is the lattice spacing. 

Physically, the choice of the time direction is related to interpreting Polyakov loop cor- 
relations as the qq potential; then the time direction is the one in which the Polyakov loop 
closes. In our study we follow Berg and Billoire to choose Polyakov loops closed in the 
iV a direction and their correlations measured along the N c direction. Then the physical 
temperature is defined to be 

T > - h (4) 

In the coming sections we shall use the physical temperature (T p ) interpretation, so we 
shall drop the subscript p, i.e. T = T p . In the last section we shall relate this to the box 
temperature (Tg) interpretation. 

After a complete study in the intermediate volume region (1 < z g < 5), Berg and Billoire 
concluded that they did not find evidence for string formation in this region, but it is 
expected to occur in larger volumes. The question this paper studies is how to understand 



the string tension measured in intermediate volumes. Is the string tension just due to finite 
size effects? Since the string tension was calculated from correlations of Polyakov loops 
closed in the N a direction, this can be considered as measuring the potential energy of a qq 
pair in a long rectangular box, with the volume V = N b 2 -N c -a 3 and at the finite temperature 
T = 1/N a a, as shown in Fig. [V]. As we know, at high temperatures, T > T c , where T c is the 
deconfining temperature, the qq is unconfined in the infinite volume limit (Nb, N c — > oo ). 
However, as the transverse size Nb is made small, i.e. r/Na > 1, one expects the side walls 
of the box would squeeze the flux lines of the qq pair to form a tube although there is no 
intrinsic string formation in this displayed in Fig. |V|. In the following we will present 

our studies about the qq flux distributions which would support the above description. 
The remaining parts of this paper are arranged as follows. Section [TJ gives the basic 



concepts of the flux measurements. Section III discusses the qq flux distributions at finite 



temperatures and finite volume effects. Section [Iy| discusses the relation between the qq 
flux distributions in a finite box and the string tension in intermediate and large volumes. 
Finally, Section [V] gives the summary. 



II. MEASUREMENTS OF THE FLUX DISTRIBUTIONS 

In our study of the SU(2) lattice gauge theory (LGT) we use the standard Wilson action, 

%)=/?E(1-Vp), (5) 
p z 

where Up is the product of link variables around a plaquette, i.e. Up{n) = U^ir^Uyin + 
//)C/7 1 (ri + i>)U~ l {n). We can measure the flux distributions of a qq pair by calculating the 
quantity flTj, 



< P(0)Pt( r -) D ^(f) > 



(6) 



< P(0)Pt(f) > 

where P{r) = |Tr Ilf^i U a (f, r), is the Polyakov loop closed in the N a direction and n^ u = 
^Tr(Up) is the plaquette variable with the orientation (fi, u), which has 6 different values, 
Gu,z/) = (2,3), (1,3), (1,2), (1,4), (2,4), (3,4). 



To reduce the fluctuations of the quantity P(0)Pt(r)D, in practical calculations we mea- 
sure the quantity M 



< P(0)Pt(f) D ^(f) > _ < P(0)Pt(f)n^(f R ) > 



(7) 



< P(0)Pt(r) > 

as the flux distribution instead of Eq. (|^), where the reference point xr was chosen to be 
far from the qq sources. This replacement does not change the measured average due to 
the cluster decomposition theorem. The six components of /' in Eq. (|^) correspond to the 
components of the chromo-electric and chromo- magnetic fields {£, B) in Minkowski space, 
i.e. 

f» v ^\{-Bl-Bl-Bl£l£l£l). (8) 
We then define the total electric energy density to be 

Pei = \[Sl + Sl + SH (9) 
and the total magnetic energy density 

p m a = \[Bl + Bl + Bl}. (10) 

In the following we will concentrate on studying the total energy and action densities pe 
and pa, which are the combinations of p e i and p ma , The total energy density is 

PE = Pel + Pma, (H) 

The total action density is 

PA = Pel - Pma- (12) 

In our measurements we transformed our flux data from lattice units to physical units by 
using the scaling relation given by Table |, which were obtained from a similar table in ref. 
H and we interpolated points in the region 2.22 < (3 < 2.50. The detail of the interpolation 
process is presented elsewhere 0. 
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III. FLUX DISTRIBUTIONS AT FINITE TEMPERATURES 



It is well known that SU(N) gauge theory has a deconfining phase transition at some 
finite temperature, T c [JITJ. This finite temperature phase transition has been studied exten- 
tively in SU(N) LGT. The transition temperature T c can be determined from Monte Carlo 



calculations, for example, the critical temperature for SU(2) was determined to be [|TT 



P c = 2.2985 ± 0.0006 (for iV a = 4). (13) 

The relation between (3 C and T c is, T c = 1/N a a((3 c ), as given by Eq. (f|). The lattice spacing 
a is a function of the coupling constant f3, such as that of Table |[ The phase transition is 
expected to occur in the infinite volume limit, i.e. Nb, N c — > oo. In this limit, for (3 < /3 C , 
the system is in the confined phase, otherwise, it is in the unconfined phase. 

We measured the qq flux distributions on lattices N a • N% ■ N c with N a = 4, N b = 5, 7, 
9 and 11, and N c = 65. Since N c 3> A^ a , A^^, as Nb/N a gets large, we expect the system 
to approach the infinite volume limit and one should see the two phases. However, for Nt, 
small, that is, Nb satisfying the condition, 

N b ~ N a , and N b a/r < 1, (14) 

the qq pair is in a finite box and one expects finite volume effects to be large. In the 
following we will study the flux distributions with (3 > /3 C and (3 < f3 c in various spacial 
volumes respectively. 



A. Flux Distributions With (3 > f3 c 

In this case the qq system approaches the unconfined phase as the volume becomes large 
(Nb — ► oo). In this phase we expect that there is no string formation. However, for small 
volume one expects finite volume effects to be large. 

In Fig. |V| we show some typical results of the energy density pe distributions in the 
region of /3 > (3 C ((3 — 2.40). They are the flux distributions on the transverse plane midway 

6 



between the qq sources with fixed separation r = 4a, and they were measured on lattices of 
4 different spatial size iV 6 =5, 7, 9, 11. For iV& small (i.e. iV 6 =5, 7) the transverse plane is 
the whole lattice. For large Nb (Nb=9, 11) the data is truncated on the margins. The signal 
is lost in the noise beyond the region shown. 

From Fig. [V] we can compare the flux distributions in boxes of different transverse size 
Nb- One can clearly see for Nb small the flux density pE at the edges of each plane have large 
values, as shown in Fig. [V| (a), (b). This implies that finite volume effects are significant 
when Nb is small because in small volumes flux lines of the qq sources would be squeezed 
dramatically by side walls of the box, so the values of pE at edges are large in this case. As 
we increase the transverse size, A^, finite volume effects become smaller, the flux density pe 
at edges decrease rapidly to zero, as shown in Fig. |V] (c), (d). 

Here we want to emphasize that the values of pe at the edges of each plane are not due 
to the reference point in Eq. fl7|). In our flux measurement we choose the reference point xr 
far from the qq sources, and we find the reference value < P{9)P\r)^^ u {xR) > in Eq. (|7|) 
is consistent with the product < P(0)P\f) >< > within errors. 

In Table |TJ we list the typical values of flux densities pe, Pa at edges of each plane in 
Fig. [V] and their corresponding errors. From this table one can explicitely see that the values 
Pe and pa at edges decrease rapidly with Nb, as we observed from Fig. |V|. However, we also 
notice that even in cases of large Nb (i.e. Nb=9, 11) the values of pa shown are non-zero 
within errors. This may be caused by a number of factors. The edge of the plane is not 
the boundary of the lattice for cases of Aj,=9, 11, the flux values may in fact vanish on the 
boundary of the lattice, but we did not calculate them there for practical reasons. Also Nb 
is perhaps not large enough (even for Nb = 11), so there are still some small finite volume 
effects. Finally our error bars were only calculated from statistical error, the actual errors 
may be larger due to systematic errors. 

To see the behaviour of the flux distributions changing with the qq separation r, in Fig. [V] 
we show the pA distribution on the transverse plane for 4 different separations, r = 3a, 4a, 
5a and 6a. The flux data were measured on the lattice 4 • 5 2 • 65 with (3 = 2.40. From this 
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figure one can see the peak values of pa decrease rapidly with the increase of r. At large 
r (i.e. r = 6a) the peak of pA almost vanishes, the flux density on the plane approachs a 
uniform distribution. 

We then calculated the center slice energy oe and action a a from our flux data, which 
are the energy and action stored in the transverse slice of unit thickness midway between 
the qq pair. The results are shown in Fig. [V[ In this figure we plot the behaviours of <je vs. 
r and a a vs. r respectively for (3 = 2.36. This shows that for A& small ge and a a do not 
decrease to zero as r increaces. However, for iV& large (e.g. A& = 9, 11), <je and a a decrease 
rapidly with r and become very small at large qq separations (i.e. r = 6a). 



B. Analysis Of Finite Volume Effects For Cases With (3 > f3 c 

To see how the finite volume effects influence qq flux distributions in a finite box, let us 
consider an electrostatic charge pair +e, — e enclosed in a similar long rectangular box as 
that of Fig. [V[ The interaction between charges is the Coulomb interaction V(r) ~ 1/r. As 
the charge separation r becomes very large, i.e. r/N^a — > oo, one can assume the electric 
field S on the middle transverse plane is uniform, and can be written as 

e = j£p (15) 

where $ is electric flux through the transverse plane, which is a constant for the Coulomb 
interaction, and (A^a) is the transverse size of the box. So the total electric field energy on 
the transverse plane is 

This shows that at large charge separations the center slice energy {(Je)c decrease with the 
behaviour (A^a) -2 as the transverse size Nf, increases, and {<Je)c vanishes as — > oo, where 
the label £ C denotes Coulomb interaction. 

Now let us return to the qq pair in the box with (3 > (3 C . As we discussed in the above 
section, the system approaches the unconfined phase for large A&. In this phase the qq 
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interaction is a screened Coulomb interaction V(r) ~ e mr jr with m the screening mass 
1~2|| . Then if the qq pair is put in a box such as that in Fig. [V], the flux $ through the 



transverse plane is not a constant, but would decrease with r and Nj,. So we expect that 
in this case the center slice energy oe for a qq pair would decay faster with Nf, than (<7e)c, 
which has the inverse square behaviour, (<Je)c ~ (N b a)~ 2 , for r — > oo. In Fig. |V|we plot our 
erg data versus the transverse size Nj, for large qq separation r. Our data are compared with 
the Coulomb behaviour (iV^a) -2 and for reference purpose the inverse quartic behaviour 
(Nbd)~ 4 . The data were measured on lattices 4 • N£ ■ 65 with iV6 =5, 7, 9, 11 and (3 = 2.40. 

From Fig. [V| one can see that for large r our oe data appears to decay faster with Nb than 
the Coulomb behaviour (iV^a) -2 as expected. This shows us that the qq interaction in the 
unconfined phase at least contains a term that decays faster than the Coulomb interaction, 
such as the screened Coulomb interaction, although the data is not good enough to determine 
the screening mass. 

In conclusion, our flux data in the region of (3 > /3 C shows that in the unconfined phase 
there is no string formation. For small transverse size iV&, the finite volume effects are large, 
the flux lines between a qq would be squeezed by side walls of the box significantly. This 
would result in a finite string tension. 



C. Flux Distributions With (3 < (3 C 

For j3 < /3 C the qq system approaches the confined phase as the volume becomes large 
(Nb — > oo). We expect string formation would occur in this phase. To see this we need to 
study behaviours of the flux distributions as a function of the qq separation r. 

In Fig. |V| we show the action density pa distribution changing with the qq separation 
r. The flux data were measured on the lattice 4 ■ ll 2 ■ 65 with (3 = 2.25. Since Nb is large 
(iVk=ll) we deleted the margins of the transverse cross section of the lattice, as we did in 
Fig. [V|. This figure should be compared with Fig. |V|, which is for (3 > (3 C . One can see Fig. [V] 
shows significantly different behaviour from Fig. 0. In Fig. the peak values of pA on the 
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plane approach a finite value as r becomes large. Even at large r (i.e. r = 6a) this peak 
still exists. However, in Fig. [V] the peak of the pa distribution almost disappears at r = 6a. 
So the flux distribution in Fig. |V| implies that intrinsic string formation occurs in the region 
of f3 < f3 c . This is not due to finite volume effects because these effects are small at large 
volumes (i.e. N b =9, 11), as we discussed in Fig. |V| and Table |TJ. 

We also calculated the center slice energy and action ce, a a from our flux data in the 
region (3 < (3 C . If string formation occurs in the confined phase, for Nj, large both oe and 
a a should approach to some finite non-zero constants when r — > oo. In Figs. [V|we plot the 
behaviours of oe vs. r and a a vs. r respectively for j3 = 2.28. The data were measured 
on lattices of various spatial size N b =5, 7, 9, 11. From this figure one can see that in all 
cases oe and a a do not decrease with r. For each fixed N b the values of oe and a a are 
almost constant as r increases. Further, for large iV& (N b = 9 ,11), where finite volume 
effects are small, both <je and a a keep as finite non-zero constants as r becomes large. This 
behaviour is totally different from that of Fig. [V], which is in the region of f3 > (3 C . So Fig. [V] 
also implies that string formation occurs. In this figure we also notice that fluctuations of 
the data are large compared to the unconfined data. This is a typical behaviour because 
confinement corresponds to disorder. 

In conclusion, our flux data in the region of (3 < f3 c provide evidence for intrinsic string 
formation in the confined phase. This string formation is not due to finite volume effects 
because these effects are small as the volume becomes large (i.e. N b = 9, 11). 

IV. SU(2) GAUGE THEORY IN FINITE VOLUMES AT T B 

As we discussed in section |, we have two ways to interpret the LGT calculations. One 
way is to identify the shortest extent N a of lattices to be the temporal size. This is a 
convenient way to study the LGT system. In previous sections we have discussed our 
calculation results in this way. However, to compare with analytical results in finite volumes 
at zero temperatures, one can use another way to interpret LGT results, that is, the longest 
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extent N c of lattices is chosen as the time direction, so that the temperature is as low as 
possible. In this case Polyakov loops closed in N a (or iV&) direction are no longer viewed 
as quark sources, they are considered to be spatial operators. Quantities such as, glueball 
mass and string tension, can be calculated from Polyakov loop correlations along the time 
direction (N c ). No matter which way we choose, quantities calculated in LGT are the same, 
the mathematics of the two ways are equivalent. In this section we shall look at the system 
in terms of this second interpretation, and reinterpret the results of previous sections. 

If we choose N c as the time direction, we are studying SU(2) gauge theory in the volume 
U = N a - N% ■ a 3 at near zero temperature Tb = 1/N c a because N c is large. For convenience, 
instead of z g defined in Eq. ([!]) we shall use another parameter, z K = y/nL = y/nN^a, to 
characterize the physical size of the volume ||, where k is the string tension measured by 
two Polyakov loops closed in N a with correlations measured along N c . This corresponds to 
Kt of ref. ||. J. Koller and P. van Baal used z g in their analytical calculations ||, so that 
they could easily compare with Monte Carlo results. B.A. Berg and A.H. Billoire used both 
z g and z K parameters and other z parameters in analyses of their LGT Monte Carlo data 
for convenience [||]. The parameter z K is equivalent to z g , For example, from the data of ref. 
one can see, z g « 1 corresponds to z K ps 0.24, and z g ~ 5 corresponds to z K ~ 1.3. In 



Table |TJ we show the correspondence of z g and z K . 



A. Results In Intermediate Volumes (0.24 < z K < 1.3) 

In this part we will show that our lattices with (3 > (3 C belong to the intermediate volume 
region. Since in this region our string tension data have large error bars, we just simply use 
the string tension data of ref. to show that our lattices with (3 > f3 c satisfy 0.24 < z K < 1.3. 
The data of ref. || are given in Table [TV], which have high statistical accuracies. 

From Table [IV] we can see the values of string tension y/Ha and the parameter z K decrease 
with the increase of Nf, and (3. Our data were measured on lattices of the size 4 • N% ■ 65 



with Nt =5, 7, 9, 11 and (3 = 2.36 and 2.40, which are similar to the lattices in Table [TV 
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We can estimate our string tension and values of z K from Table |IV|. For example, the 
estimated string tension for the lattice 4 ■ 5 2 • 65 with (3 = 2.36 could be interpolated from 
the string tension measured on lattices 4 ■ N£ ■ 64 with Nb=A, 6 and (3 = 2.36. From 
Table [TV] one can see that for A^=5 the string tension satisfies, 0.1593 < y^a < 0.2475, and 
0.80 < z K = \J~KaNb < 1.24, which is in the intermediate volume region, 0.24 < z K < 1.3, as 



shown in Table |TT[ By comparing our lattices with Table [TV| in this way, we find that our 
lattices with (3 = 2.36 and 2.40 for Nb =5, 7, 9, 11 are all in the intermediate volume region, 
with the "box temperature" approaching zero, i.e. Tb — > 0. 

As we discussed in section pit for (3 > (3 C there is no intrinsic string formation. However, 
since finite volume effects are large when the transverse size iV& is small, this results in the 
observed string tension in these cases. As Nj, becomes larger, finite volume effects becomes 



smaller, the string tension becomes smaller. This is confirmed in Table [Ty|. Since the lattices 
with (3 > (3 c 'm our study all belong to the intermediate volume region, in these cases we find 
there is no intrinsic string formation, and the string tension is due to finite volume effects. 
In general, we expect that the results apply to the whole intermediate volume region in the 
zero temperature limit (Tb — > 0). 



B. Results In Large Volumes (z K > 1.3) 

For (3 < (3 C we can easily extract the string tension, which is given in Table |V|. In this 
table we also show the values of z K for lattices with f3 < j3 c . We can see that most lattices 
listed in Table [V] satisfy z K > 1.3, with one case in the critical region, z K ~ 1.3 for (3 = 2.28 
and Nb = 5. From Table [TTT| we know that lattices satisfying the condition, z K > 1.3, belong 
to the large volume region in the zero temperature limit (Tg — > 0). So our lattices with 
[3 < j3 c , as listed in Table [V], are in this region. 

As we discussed in section [111 C[ a system with /3 < (3 C is in the confined phase in the 
infinite volume limit. In this case we have shown that intrinsic string formation occurs. As 
the transverse size Nb is made small (e.g. iV& = 5), we also observed some finite volume 
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effects, however, it does not change the nature of qq confinement. We have shown that 
the lattices with f3 < j3 c in our study all belong to the large volume region in the zero 
temperature limit (Tg — > 0). From our study we have found evidence for intrinsic string 
formation in these cases, which is not due to finite volume effects. We expect that the 
intrinsic string formation occurs in the whole large volume region. 



We used two ways to interpret the LGT results calculated on lattices N a • Nj* ■ N c with 
the geometry iV a < Nf, <C N c . First, we chose the convenient way to interpret our data, 
by identifying iV a to be the time direction we can study the system at finite temperatures 
(T p = 1/N a a) in the volume V = N£ ■ N c ■ a 3 . We analysed the flux distributions of a qq 
pair in terms of this interpretation. If we choose N c to be the time direction we can study 
the system at near zero temperature (Te = 1/N c a) in the volume U = N a ■ N£ ■ a 3 . We find 
that for (3 > (3 C there is no intrinsic string formation. Our lattices with (3 > f3 c belong to 
the intermediate volume region. We then established that the origin of the string tension 
measured by Berg and Billoire [g] in this region is related to finite sizes of the lattice. For 
(3 < j3 c we find clear signals for intrinsic string formation. Our lattices with (3 < (3 C are in 
the large volume region, but near the borderline of intermediate volumes and large volumes. 
This is just beyond the volume region investigated by Berg and Billoire. 

Another remark on string formation is that for higher units of 't Hooft electric flux one 
expects that string formation predicts the relation [|], 



V. SUMMARY 



E n = V^E_ 



(n = l,2,...), 



(17) 



where E n is the energy of the n unit flux. However, both analytical results Q and LGT 
calculations produce a different behaviour in the intermediate volume region, 



E, 



"n 



nE Xl 



(18) 
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where n = 1, 2, 3 for SU(2). Even in the large volume region some Monte Carlo data of ref. 
still show the behaviour of Eq. (|18|). We speculate that Eq. (|I7|) may be due to lattice 
artifacts since the origin of v^2 and v^3 in this equation are for the planar diagonal and 
volume diagonal on lattices. 

Our study supplement the global study of Berg and Billoire ||, and provide evidence 
that string formation does not occur in intermediate volumes, but occurs in large volumes. 
To verify Eq. fll7D or Eq. (|18|) we need to measure the energy E n of higher unit flux in the 
large volume region. 
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FIGURES 



The qq colour sources in a long rectangular box. The flux lines between the qq pair are 
confined by side walls and have a flux tube form. 

The energy density p E flux distributions in the region of (3 > (3 C on the plane midway 
between the qq sources of separation r = 4a. The flux data were measured on lattices 
4 • iV 2 ■ 65 with (3 = 2.40 for various spatial sizes, (a) N b = 5, (b) N b = 7, (c) N b = 9, and 
(d) N b = 11. These flux data are measured in the physical unit Gev/ fm 3 . 

The action density pa flux distributions in the region (3 > f3 c on the plane midway 
between the qq sources of various separations (a) r = 3a, (b) r = 4a, (c) r = 5a and (d) 
r = 6a. The flux data were measured on lattices 4 • 5 2 ■ 65 with (3 = 2.40, The data are 
measured in the physical unit Gev/ fm 3 . 

The behaviour of center slice energy a E and action a a verses the qq separaction r in the 
region of (3 > j3 c , (a) a E vs. r, (b) a a vs. r. The data were measured on lattices 4 • • 65 
with P = 2.36 for various spatial size, N b = 5 (circles), N b = 7 (squares), N b = 9 (triangles), 
N b — 11 (diamonds). The data are in the physical unit Gev/ fm. 

The plot of ge in the region of (3 > (3 C verses the transverse size of lattices N b at large qq 
separation r = 6a. The solid line is the Coulomb behaviour (N b a)~ 2 , the dashed line is the 
inverse quartic behaviour (N b a)~ 4 , both are normalized to the data at N b = 5. The data 
were measured on lattices 4 • iV b 2 • 65 with N b =5, 7, 9, 11. and (3 = 2.40. 

The action density pa flux distributions in the region of (3 < (3 C on the plane midway 
between the qq sources of various separations (a) r = 3a, (b) r = 4a, (c) r = 5a and (d) 
r = 6a. The flux data were measured on lattices 4 • ll 2 • 65 with j3 = 2.25, The data are 
measured in the physical unit Gev / fm 3 . 
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The behaviour of center slice energy oe and action a a verses the qq separaction r in the 
region of f3 < (3 C) (a) a E vs. r, (b) a A vs. r. The data were measured on lattices 4 • iV b 2 • 65 
with f3 = 2.28 for various spatial size, iVj, = 5 (circles), Nb = 7 (squares), Nb = 9 (triangles), 
iVf, = 11 (diamonds). The data are in the physical unit Gev/fm. 
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TABLES 

TABLE I. The correspondence of the lattice spacing a and the coupling constant j3 for SU (2) 
LGT. 

P 2.25 2.28 2.36 2.40 

a (fm) 0.1966 0.1785 0.1374 0.1203 



TABLE II. The typical measured flux density values (errors) of pe and pa in the region of 
(3 > j3 c on the edges of the transverse plane midway between the qq, as shown in Fig. The flux 
data are in the physical unit Gev / frrfi . 



4 • Nj; • 65 


r = 4a 


N h 


5 


7 


9 


11 


P = 2.40 


Pe 


2.5 (5) 


0.86 (24) 


0.11 (18) 


0.31 (17) 




PA 


17.0 (6) 


2.9 (3) 


1.5 (2) 


1.0 (2) 



TABLE III. The correspondence of the two parameters z g = m(0 + )Nf,a and z K = y/nN^a, 
which are obtained from data in ref. ||. 



small volume 


intermediate volume 


large volume 


Zg < 1 


1 < Zg < 5 


Zg>5 


z K < 0.24 


0.24 < z K < 1.3 


z K > 1.3 
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TABLE IV. The values of the strin 


g tension data and z K 


on lattices of size 4 • 


N% ■ 64 with 




=4, 6, 8 and (3 = 2.36, 2.38 and 2.41. 


. The data are quoted from ref. ||. 




4 • N(f ■ 64 




N b 


4 


6 


8 


P = 


- 2.36 \fka 


0.2475 (18) 


0.1593 (39) 


0.1091 (39) 




z K = N h ^J~Ka 


0.99 (1) 


0.96 (2) 


0.87 (3) 


P = 


= 2.38 yJ~Ka 


0.2424 (23) 


0.1357 (44) 


0.0732 (81) 




Z K = Nhx/KCL 

fx, U V 


0.97 (1) 


0.81 (4) 


0.59 (6) 


P = 


= 2.41 ^/ko 


0.2264 (14) 


0.1314 (51) 








0.91 (1) 


0.79 (3) 






TABLE V. The string tension data 


k and the values of the 


z K parameter on lattices of the size 


4- 


N£ ■ 65 with N b =5, 7, 9, 11 and f3 = 


2.25 and 2.28. 






4 • iV 6 2 • 65 


N b 


P 






z K = N b y/na 


5 


2.25 


0.300 (22) 




1.5 (1) 




2.28 


0.251 (34) 




1.3 (2) 


7 


2.25 


0.288 (31) 




2.2 (2) 




2.28 


0.258 (31) 




1.8 (2) 


9 


2.25 


0.262 (84) 




2.4 (8) 




2.28 


0.226 (67) 




2.0 (6) 


11 


2.25 


0.307 (48) 




3.4 (5) 




2.28 


0.232 (64) 




2.6 (7) 
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